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We show that, on the contrary to the usual claims, fully supersymmetric CP violations in the 
kaon system are possible through the gluino mediated flavor changing interactions. Both ex and 
Re(e'/e_ft- ) can be accommodated for relatively large tan/3 without any fine tunings or contradictions 
to the FCNC and EDM constraints. 
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Until this year, the only CP violation observed was 
in Kl — > 27r [Q, which could be attributed to AS = 2 
K° — K° mixing. The mixing parameter ek is accurately 
measured by now : e K = e"/ 4 (2.280 ± 0.013) x 10" 3 ||. 
Recent observation of Re(e'/e#r) by KTeV collaboration, 
Eje(e'/e K ) = (28±4) x 10" 4 ||, nicely confirms the earlier 
NA31 experiment (§ Re{e'/e K ) = (23 ± 7) x 10" 4 . This 
nonvanishing number indicates unambiguously the exis- 
tence of CP violation in the decay amplitude (AS = 1). 
These two parameters quantifying CP violations in the 
kaon system can be accommodated by the KM phase 
in the Glashow-Salam- Weinberg's standard model (SM). 
The SM prediction for the latter is about 5 x 10 -4 and 
lies in the lower side of the data, although theoretical 
uncertainties from nonperturbative matrix elements and 
the strange quark mass are rather large [pi. 

However, it would be interesting to consider a possi- 
bility that these CP violations have their origin entirely 
different from the KM phase in the SM, in particular in 
the framework of various extensions of the SM including 
supersymmetric models f|. In the minimal supersym- 
metric standard model (MSSM) considered in this work, 
there are many new CP violating phases that fall into 
two categories : phases with flavor preserving (FP) and 
flavor changing (FC), each of which is constrained by 
electron/neutron electric dipole moments (EDM's) and 
the €k , respectively. Recently, it was shown that the FP 
and CP violating phases in A t and /i in the more min- 
imal SUSY model do not generate enough ck @ @] or 
new phase shift in B° - B° (?J ||, although they can 
lead to a large direct CP asymmetry in B —* X s j upto 
~ ±16% if charginos and stops are light enough ||. How- 
ever, another class of CP violating phases in the flavor 
changing quark-squark-gluino vertices in the MSSM may 
be relevant to CP violations in the K meson system. 

In this letter, we show that all the observed CP violat- 
ing phenomena in the kaon system in fact can be accom- 
modated in terms of a single complex number (8f 2 )LL 
that parameterizes the squark mass mixings in the chi- 
rality and flavor spaces for relatively large tan /3 without 
any fine tuning or any contradictions with experimental 
data on FCNC, even if 5km = 0. We assume that CKM 
matrix is real in the most of this letter for simplicity and 



maximizing the effect of our mechanism. The case where 
the KM phase is nonzero is discussed in brief, and the 
details will be given elsewhere Jl(J . 

In order to study the gluino (photino) mediated flavor 
changing phenomena in the quark (lepton) sector such as 
ArriKi^K and Re^'/e^) or lepton flavor violations, it is 
convenient to use the so-called mass insertion approxima- 
tion (MIA) |ll[] . The quark-squark-gluino vertex is flavor 
diagonal in the MIA, and the flavor/chirality mixing oc- 
cur through the insertion of (<5^)ab, where i,j = 1,2,3 
and A, B = L, R denote the flavors of the squark under 
consideration and the chiralities of its superpartner. The 
superscript denotes that the down type squark mass ma- 
trix is involved. The parameters (Sfj) ab characterize the 
size of the gluino-mediated flavor changing amplitudes, 
and they may be CP violating complex numbers, in gen- 
eral. In the following, diagrams involving charged Higgs, 
chargino and neutralino will be ignored, since they are 
suppressed by a w /a s compare to the gluino-squark loops 
unless gluino/squarks are very heavy. This should be a 
good starting point for studying the SUSY FCNC/CP 
problems. 

Now, if one saturates Atuk and e# with (Sf 2 )LL alone, 
the resulting Re(e /crt) is too small by more than an 
order of magnitude, unless one invokes some finetuning 
|l2f . Recently, Masiero and Murayama showed that this 
conclusion can be evaded in generalized SUSY models 
|fl3f with a few reasonable assumptions on the size of 
the (#12) lr and the relations between the down quark 
Yukawa couplings and the CKM mixings. But they did 
not consider possibility to generate ck from a supersym- 
metric phase in Ref . [1J] , and also predict too large neu- 
tron EDM which is very close to the current upper limit. 

In the following, we show that there is another generic 
way to saturate e /ck in supersymmetric models if 
|/itan/?| is relatively large, say ~ 10 — 20 TeV. More- 
over, both ck and Re(e /ex) can be generated by a sin- 
gle CP violating complex parameter in the MSSM. In 
other words, fully supersymmetric CP violations are pos- 
sible in the kaon system. The argument goes as follows 
: if l(<^i 2 )iil ~ O(10- 3 - 10- 2 ) with the phase - 0(1) 
saturates ck , this same parameter can lead to a sizable 
Re(e'/e/f) through the (Sf 2 )LL insertion followed by the 



FP (LR) mass insertion, which is proportional to 

(8 d 22 ) LR = m s (A* - M tan/?)/m 2 ~ 0(HT 2 ), 

where rh denotes the common squark mass in the MIA. 
It should be emphasized that the induced (Sf 2 y^ R = 
{&\2)ll x {&22)lr is different from the conventional 
(^12) lr in the literature. The loop functions for these 
two LR insertions are different with each other in general. 
The LR mixing (Su) 1 ^ induced by (<5f 2 )z,L is typically 
very small in size ~ O(10 -5 ), but this is enough to gener- 
ate the full size of Re(e'/e#-) as shown below. Our spirit 
to generate supersymmetric Re(e /^k) is different from 
Rcf . [^3) , where the LR mass matrix form is assumed to 
be similar to the Yukawa matrix so that they predict the 
neutron EDM to be close to the current upper limit. On 
the other hand, our model does not suffer from the EDM 
constraint at all, as shown below. 

Let us first consider the gluino-squark contributions to 
the K° — K° mixing due to two insertions of (Sf 2 ) ll- The 
corresponding AS* = 2 effective Hamiltonian is given by 

H cfi (AS = 2) = C~d a Llti s a L 2£V*£ 

with the Wilson coefficient C\ being 
,2 
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Here, x — m~/rh 2 and the loop functions fe(x) and fe(x) 
are given in WM . The double mass insertion diagrams of 
FC LL followed by the FP LR generate another oper- 
ators : Q 2 — d^s R d^s R and Q3 = d^s R d^s R , whose 

Wilson coefficients are proportional to [(612)^] 2 . Since 
their effects on ArtiK and ejf are negligible, we do not 
show them here explicitly although their effects have been 
included in the numerical analyses. 

Now we turn to the AS = 1 effective Hamiltonian 
Ti. e s(AS = 1) = X^=3 CiOi. The sdg operator O s which 
is relevant to Re(e /ex) is defined as 
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and other four quark operators Oi=3,...,6 an d the corre- 
sponding Wilson coefficients from C3 to C$ with a single 
mass insertion are available in the literature E3 . One 
has to remind that C^... t a's are proportional to (Sf 2 )LL, 
whereas C% is given by a linear combination of (Sf 2 )LL 
and (S 12) lr- This (S^ 2 )lr dependent part in C$ is pro- 
portional to m g jm s , and thus is very important for gen- 
erating Re(e /ek) even if (Sf 2 )LR is fairly small. 

If we consider the penguin diagram Fig. [l] with the 
double mass insertion, the Wilson coefficient Cs is given 
by 



a 



(2 ) _ a s nig d ind 

- ^2 \°12)LR 

m, m s 



C 1 m[ 2) (x)+C 2 M {2) (x) , (3) 



m[ 2 \x) = 



Mi 2 \x) = 



3/2, C 2 
2(3 - 



= — 1 / 6 and 
3x 2 + (1 + 4a; 



i log a;) 



(x- l) 5 

9x — 9x 2 — x 3 + 6x(l + x) logx 



1 



(x-iy 



(4) 



The contributions of photon penguin and Z penguin di- 
agrams with the double mass insertion are negligible as 
in the case of the single mass insertion jl3 . 

Now we are ready to calculate the SUSY contributions 
to Avuk^k and e'/ejc using the AS = 1,2 effective 
Hamiltonians obtained above and the following expres- 
sions Hal: 



Atok(SUSY) = 2ReMi2, 
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where 2mifM 1 * 2 = 

the isospin amplitudes defined as 
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({irir) r \K ). In the numerical analysis, we use 
the same parameters as in Ref. |l5| with m s (2GeV) = 
130MeV. The corresponding SM prediction for 
Re(e /ex) — 5.7 x 10~ 4 . We vary the modulus and 
the phase of (5f 2 )LL as indepent parameters, and select 
those points which satisfy Am^(SUSY) < Amx(exp) 
and \ck (SUSY) — e^- (exp)| < la. Then, for these points, 
we plot e' /ex in Figs. 2 (a)-(d) as functions of the mod- 
ulus r [(a) and (c)] and the phase ip [(b) and (d)] of the 
parameter (Sf 2 )LL = re %ip for the common squark mass 
rh = 500 GeV. The upper (lower) rows correspond to 
A~ s = (A s - fi* tan/?) = -10 (-20) TeV. Different x's (= 
0.3, 1.0, 2.0) are represented by the solid, the dashed 
and the dotted curves, respectively. If we choose the 
opposite sign for A s , the phase of the relevant ($i 2 )ll 
should be shifted by 180° in order that we have correct 
sign for ejf. From Figs. 2 (a) and (b), it is clear that both 
ex and Re(e'/e#-) can be nicely accommodated with a 
single complex number (Sf 2 )LL with ~ 0(1) phase in our 
model without any difficulty, if \fi\ and tan/3 is relatively 
large so that \A S \ becomes a few tens of TeV. If the com- 
mon squark mass rh differs from 500 GeV, the A t should 
be multiplied by (rh in GeV/500) 2 for the fixed x. 

Let us consider the neutron EDM constraint. The 
FP LR mass insertion in the gluino-squark diagram con- 
tributes to the neutron EDM. The effective Hamiltonian 
for the neutron EDM is given by @ H e //(EDM) = 
Ei=i Cf ra O l , where the CVs are defined as 
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In our model, the FP LR mass insertion in the gluino- 
squark diagram contributes to the neutron EDM with 
the following Wilson coefficients : 

Cf- = -l e -^Q d §lm(Sf 1 ) LR 

where (Sf^LR = rnd(A d — /itan/3)/m 2 , and 
(1) l + 4x-5x 2 + 2x (2 + x) log(x) 

5 (x) = 2M + ,) 4 ' (8) 

m , . 2 (-ll + lOx + x 2 ) - (9 + 16x-a; 2 ) log(x) 

; (x) = 2 . 

3(-l + x) 4 

and Cf dm = 0. Our expression for Cf dm confirms the re- 
sult obtained in Ref. H] , and the result for C| dm is new. 
The renormalization group (RG) running effect and the 
final formula for the neutron EDM can be found in jl6| . 
We found that the EDM constraint is very strong. If 
the universality of the trilinear couplings are assumed 
(Ad = A s ), then LR should be essentially real. Still 
a single CP violating phase in (5f 2 )LL can generate right 
amounts of ck and Re(e'/ e if) without any fine tuning. 
Even if we relax the universality condition Ad = A s , the 
result is basically the same, since we are in the regime of 
large /i tan f3 and its phase is constrained by the neutron 
EDM irrespective of Ad = A s as long as |A*,sl ^ 1 TeV. 
Note that we are not assuming any specific flavor struc- 
tures in the A terms at all, unlike many other models in 
Ref. §. 

It should be worthwhile to emphasize the importance 
of the FP (LR) mixing (8 22 )lr in our study. If we ignored 
this effect and considered the flavor changing (LL) and 
(RR) mixings simultaneously for example, we could get 
both 6k and Re(e' /ck), but significant amounts of fine 
tunings are unavoidable. The ratio of the magnitues of 
(Sf 2 )LL and (Si 2 )rr should be 0(1O~ 3 ) in order that we 
explain the large experimental data for Re(e' /ex)- Then, 
the contributions from (LL) 2 [ or (RR) 2 , whichever the 
larger one ] and (LL) x (RR) terms should cancel with 
each other within a part in 10 3 in order to reproduce the 
experimental value for 6k , thus requiring substantial fine 
tuning @]. 

If the simplifying assumption of the real CKM matrix 
is relaxed, there will be additional constributions to €k 
and Re(e /ex) from the SM and other SUSY loop dia- 
grams. If we assume that only the gluino-squark contri- 
bution considered above is comparable with the SM con- 
tribution, it would be possible that the €k is mainly dom- 
inated by the KM phase contributions, but the Re(e I^k) 
has significant contributions from the induced (S^ 2 )lr as 
disucssed in this letter ]Io| . 

The best discriminant between our model and the 
SM model would be probably the branching ratios for 
K — ► Tivv and CP violations in B decays. The branch- 
ing ratio for the decay Kl ~ > ir°i/D is essentially zero 



in our model, since it is purely CP violating but there 
is no appreciable CP violation in the s —* dvv ampli- 
tude through gluino loop diagram |Q. On the contrary, 
K + — > n + vv involves both CP conserving and CP violat- 
ing amplitudes, and the corresponding branching ratio in 
our model ranges over (8.9±7.2) x 10" 11 , compared to the 
SM prediction : B(K+ -> n+vD) SM = (7.7±3.0) x 1CT 11 
|l7f . Our predictions can be changed by two ways, how- 
ever. The chargino-upsquark loop contributions con- 
tribute to K — > txvv through the enhanced sdZ pen- 
guin vertex, but one still expects the branching ratio for 
Kl — * n°vi> to be smaller than the SM predictions Jl8| . 
Also, if our assumption on the real CKM matrix is re- 
laxed, the KM phase will contribute to the K — > irvv, 
but the predictions will differ from the SM, since there is 
additional contribution to €k so that the CKM elements 
are less constrained (HI. If we stick to 5km = 0, the 
CP violations in B system will be very different from the 
SM predictions. However, the gluino-squark loop contri- 
butions to the B system is governed by new parameters, 
(5f 3 )AB and (S 23 )ab, which are independent of (5f 2 )AB 
we considered here. Therefore we cannot make defninte 
predictions for B decays. Generically the situation could 
be very different from the SM case (T^| . 

Now let us consider a typical size of (5f 2 )hL in the 
MSSM. The answer to this question would depend on 
the models for soft SUSY breaking and how to solve 
the SUSY flavor problem. If one invokes approximate 
(abelian or nonabelian) flavor symmetry in order to solve 
the SUSY FCNC problem, the natural size of (5f 2 ) LL 
could be order of ~ (A 3 — A 4 ) ~ 10~ 3 with order 0(1) 
phase (where A = sin^c = 0.22) pC(| . This is usually 
referred as the alignment in contrast to more popular 
universality. This is in fact the twisted version of the so- 
called SUSY FCNC and SUSY e K problem, saying that 
the gluionic SUSY contribution to €k could be generically 
very large unless the squarks are degenerate and/or the 
mass matrices of quarks and squarks are almost aligned 
in the flavor space. In this context, our result on e /eK 
is nothing but to say that SUSY 6k problem in gen- 
eral implies the SUSY e problem for relatively large 
|/itan/3| - O(10) TeV. 

In conclusion, we showed that both ck and Re(e' I^k) 
can be accommodated with a single CP violating 
and flavor changing down-squark mass matrix elements 
\(&yi)ll ~ 10~ 3 ] without any fine tuning or any conflict 
with the data on FCNC processes, if |/ztan/3| ^10 — 20 
TeV with a scale factor [?ri(in GeV)/500] 2 for the fixed 
x. Our mechanism utilizes this FC LL mass insertion 
along with the FP LR mass insertion propotional to 
(S 22 )lr ~ 10~ 2 . The latter is generically present in any 
SUSY models including the MSSM, and thus there is 
no fine tuning in our model for accommodating both 6k 
and Re(e'/ex) in terms of a single (S^ 2 )ll- It is straight- 
forward to extend our mechanism including the nonva- 
nishing KM phase, (8f 2 )hR and/or (S^ 2 )rr- One can 
also consider our mechanism in the more minimal SUSY 
model, where (S 22 ) LR is proportional to mf,(Af, — /itan/3) 
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so that At — /itan/3 may be lowered significantly. All 
these finer details including B(K — > irvv) will be dis- 
cussed elsewhere in the forthcoming publication [[l(| . 
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FIG. 1. Feynman diagram for AS = 1 process. The cross 
denotes the flavor changing (LL) and the flavor preserving 
(LR) mixings, respectively. 
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FIG. 2. Re(e jtK) as a function of the modulus r [(a) and 
(c)] and the phase <p [(b) and (d)] of the parameter (<5ll)i2 
with A s to be -10 TeV ((a),(b)) and -20 TeV ((c), (d)). 
The common squark mass is chosen to be m = 500 GeV, and 
the solid, the dashed and the dotted curves correspond to 
x — 0.3, 1.0, 2.0, respectively. 
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